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Abstract —This paper studies a three-node network in which an 
intermediate nano-transceiver, acting as a relay, is placed between 
a nano-transmitter and a nano-receiver to improve the range 
of diffnslon-hased molecular communication. Motivated by the 
relaying protocols used in traditional wireless communication sys¬ 
tems, we study amplify-and-forward (AF) relaying with fixed and 
variable amplification factor for use in molecnlar communication 
systems. To this end, we derive a closed-form expression for the 
expected end-to-end error probability. Furthermore, we derive a 
closed-form expression for the optimal amplification factor at the 
relay node for minimization of an approximation of the expected 
error probability of the network. Our analytical and simulation 
results show the potential of AF relaying to improve the overall 
performance of nano-networks. 

I. Introduction 

Molecular communication (MC) is a biocompatible ap¬ 
proach for enabling communication among so-called nano¬ 
machines by exchanging information via molecules. Integrat¬ 
ing communication capabilities expands the potential function¬ 
ality of individual nano-machines such that communities of 
them, so-called nano-networks, can execute collaborative and 
challenging tasks in a distributed manner. Sophisticated nano¬ 
networks are expected to have various biomedical, environ¬ 
mental, and industrial applications ||T]. 

Diffusion-based MC is a passive and energy-efficient ap¬ 
proach for communication among nano-machines where the 
transportation of molecules from a nano-transmitter to a nano¬ 
receiver in a fluid environment relies on free diffusion only; 
no additional infrastructure is required. However, one of the 
main drawbacks of diffusion-based MC is its limited range 
of communication, since the propagation time increases and 
the number of received molecules decreases with increasing 
distance. This makes communication over larger distances 
challenging. 

One approach from conventional wireless communications 
that could be adapted for MC to aid communication with 
distant receivers is the use of relays. In fact, the relaying of 
information is already used by nature. In particular, cascades of 
signal amplification are a common way to relay information 
from the exterior of a cell to its interior Q. For example, 
the binding of a signal molecule to a cell surface receptor 
may activate a number of G protein molecules, each of which 
activates in turn a molecule of adenylyl cyclase, resulting 
in an excessive number of cAMP molecules. Subsequently, 


each cAMP molecule activates a protein kinase, which in turn 
generates several copies of a specific enzyme. Each enzyme 
can trigger a chemical reaction, resulting in a large amount 
of enzymatic product that can ultimately change the behavior 
of the cell, cf. Fig. For example, one rhodopsin pigment 
cell (the cells of the eye responsible for interpreting light and 
dark), when excited by a photon, can trigger the release of 
10® enzymatic products, i.e., cGMP molecules, so there is an 
amplification of five orders of magnitude. 

Multihop relaying among nano-machines has been studied 
in the existing MC literature; see 0-II3- In Q and 0^ a 
diffusion-based multihop network between bacteria colonies 
was analyzed, where each node of the network was formed 
by a population of bacteria. In Q and the design and 
analysis of repeater cells in Calcium junction channels were 
investigated. In Q, the rate-delay trade-off of a three-node 
nanonetwork was analyzed for a specific messenger molecule, 
polyethylene, and network coding at the relay node. The use 
of bacteria and virus particles as information carriers in a 
multihop network was proposed in Q and Q, respectively. 
Most recently, the authors in m investigated a two-hop 
network, where a relay node is installed to maximize the 
concentration of molecules at the destination. However, most 
prior works consider the transmission of a single symbol (or 
bit) and, consequently, the effect of intersymbol interference 
(ISI), which is unavoidable if multiple symbols (bits) are 
transmitted, is not taken into account. An exception is our 
work in ini, |Tg, where the performance of a decode- 
and-forward (DF) relaying protocol is studied. However, the 
complexity associated with full decoding at the relay may 
not be affordable in certain applications. On the other hand, 
installing a simple amplifier between a nano-transmitter and a 
nano-receiver may be easier to realize. 

In this paper, we assume that the transmitter nano-machine 
emits multiple random bits and investigate two different re¬ 
laying protocols: variable-gain and fixed-gain amplify-and- 
forward (AF) relaying. In variable-gain AF relaying, the relay 
nano-machine amplifies the signal received from the nano¬ 
transmitter by a variable amplification factor, which may vary 
in each bit interval, to forward it to the nano-receiver. In fixed- 
gain AF relaying, the amplification factor at the relay nano¬ 
machine is constant for all bit intervals. The main contributions 
of this paper can be summarized as follows: 
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Fig. 1. Example for amplification of a signal at different stages of a typical 
cell-signaling process. Reproduced from j^. 


1) We derive closed-form expressions for the expected end- 
to-end error probability of variable-gain and fixed-gain 
AF relaying. 

2) We optimize the performance of variable-gain AF re¬ 
laying by deriving a closed-form expression for the 
optimal amplification factor at the relay nano-machine. 
We propose a mechanism for variable-gain AF relaying 
where the relay nano-machine estimates the optimal 
amplification factor in each bit interval. 

3) For fixed-gain AF relaying, the optimal amplification 
factor is derived by averaging the optimal amplification 
factor of variable-gain AF relaying over all bit intervals. 

The rest of this paper is organized as follows. In Section 
[m we introduce the system model and preliminaries of the 
error rate analysis. In Section we derive the expected 
error probability of the two-hop network for both AF relaying 
protocols and the optimal amplification factor. Numerical 
results are presented in Section [rV| and conclusions are drawn 
in Section IV] 

II. SYSTEM MODEL AND PRELIMINARIES 

In this section, we introduce the system model and some 
preliminaries required in the remainder of the paper. 

A. System Model 

In this paper, we use the terms “nano-machine” and “node” 
interchangeably to refer to the devices in the network, as the 
term “node” is commonly used in the relaying literature. We 
assume that a source (S') node and a destination {D) node are 
placed at locations (0,0,0) and (a:£>,0,0) of a 3-dimensional 
space, respectively. The relay (i?) node is placed in the middle 
between node S and node D along the x-axis. We assume 
that node D and node R are spherical in shape with fixed 
volumes (and radii) Vd {ro) and Vr {tr), respectively, and 
that they are passive observers such that molecules can diffuse 
through them. Eor example, small, uncharged molecules, such 
as ethanol and urea, can enter and leave a cell by passive 
diffusion across the plasma membrane, see 0. 

We assume that there are two distinct types of messenger 
molecules namely, type Ai and type A 2 , and that relay R 
can detect type Ai molecules, which are released by node S, 
and emits type A 2 molecules, which are detected by node D. 


The number of molecules released of type Af^ f € {1, 2}, is 
denoted as Naj, and the concentration of type Af molecules 
at the point defined by vector r at time t in molecule • m“^ 
is denoted by CAf{f,t)- We assume that the movements of 
individual molecules are independent. 

The information that is sent from node S to node D 
is encoded into a binary sequence of length L, Wr = 
{VFs[l], W5[2], ..., VF 5 [L]}. Here, VFs[j] is the bit transmitted 
by node S in the jth bit interval with Pr(VFs[i] = 1) = 
and Pr(Ws[j] = 0) = Pq = 1 —Pi, where Pr(-) denotes prob¬ 
ability. Eor compactness, we denote a subsequence transmitted 
by node S by = {VFs[a],..., fEs[6]}. The information 

bit detected at node D in the jth bit interval is denoted 
by Wulj]- We adopt ON/OFF keying for modulation and 
a fixed bit interval duration of T seconds. Node S releases 
Nax molecules at the beginning of the bit interval to convey 
information bit “1”, and no molecules to convey information 
bit “0”. Furthermore, we consider a full-duplex scheme, where 
reception and transmission occur simultaneously at the relay 
node, i.e., in each bit interval, relay R receives a signal 
from node S, and releases an amplified version of the signal 
received in the previous bit interval to node D. 


B. Preliminaries 


In the following, we consider a single communication link 
between a transmitting node n and a receiving node q, and 


review the reception mechanism at node q, cf. 113|. The two 
nodes can form the source-to-relay, {n = S,q = R), or relay- 
to-destination, (n — R,q — D), hop in our system model. 
For conciseness of presentation in this subsection, we drop 
the subscript / and denote the type of molecule released by 
node n and detected at node q by A. 

The independent diffusion of molecules through the envi¬ 
ronment can be described by Tick’s second law as 

dCA{r,t) 


dt 


= DAS7^CA{r,t), 


( 1 ) 


where Da is the diffusion coefficient of A molecules in —. 
Assuming that node n is an impulsive point source, and emits 
Na molecules at the point defined by vector r„ into an infinite 
environment at time t' — 0, the expected local concentration 
at the point defined by vector r and at time t is given by 114 
Eq. (4)] 


CA{r,t) = 


Na 


■ exp 


\r-rn\ 

4Da 


( 2 ) 


It is shown in that the number of molecules observed 
within the volume of node q, Vq, at time t due to one emission 
of Na molecules at r„ at t' = 0, i.e., N^’'^{t), can be 
accurately approximated as a Poisson random variable (RV) 
with time-varying mean given by 


Nil‘^2it)=CA{rq,t)Vq, 


(3) 


where Pq is the vector from the origin to the center of 
node q. Eq. 0 was derived under the uniform concentration 
assumption, i.e., it was assumed that node g is a point observer 




























or that the concentration throughout its volume is uniform and 
equal to that at its center; see GD- This is a valid assumption 
when node q is far from node n. The probability of observing 
a given A molecule, emitted by node n at t' = 0, inside Vg at 
time t, i.e., is given by (|^ when setting Na = 1 , 

i.e., 


P. 


(n,g) 


ob,A 


(t) = 




(47rZ3^f)3/^ 


exp 


4Da 


(4) 


For reception, we adopt the family of weighted sum de¬ 
tectors introduced in |T3), where the receiving node takes 
multiple samples within a single bit interval, and adds up 
the individual samples with a certain weight assigned to each 
sample. When detection is required at node g, i.e., when 
{n = R,q — D), then the sum is compared with a decision 
threshold. For simplicity, we assume equal weights for all 
samples. The decision in the jth bit interval is given by 
Eq. (37)] 

0 otherwise. 


A. Variable-Gain AF Relaying 

Node S emits type Ai molecules, which have diffusion 
coefficient Da^ and can be recognized by relay node R. 
In response to the received Ai molecules, the relay emits 
type A 2 molecules having diffusion coefficient Da 2 - Node S 
releases a fixed number of molecules, Na^ , and no molecules 
to transmit bit “ 1 ” and bit “ 0 ” at the beginning of a bit interval, 
respectively. However, the number of molecules released by 
relay node R at the beginning of a bit interval, Na 2 [•], varies 
in each bit interval. In particular, Na 2 [■] depends on the time- 
varying and random number of Ai molecules observed at the 
relay node in the previous bit interval. 

Assuming full-duplex relaying at the relay node, the trans¬ 
mission of information bit W 5 [j] from node S to node D 
has two phases. In the first phase, at the beginning of the 
jth bit interval, node S transmits information bit IFs[j], and 
node R releases Na 2 [j] molecules concurrently to amplify and 
forward the message received in the previous bit interval, i.e.. 




( 8 ) 


where is the detection threshold of node q, and we assume 
that node q takes M equally-spaced samples in each bit 
interval. The sampling time of the mth sample in the jth 
bit interval is t{j, m) = (j — 1)Tb + tm, where tm = mto 
and fo is the time between two successive samples. We 
assume that Iq is large enough for samples to be independent. 

{t{j,m)) is a Poisson RV with mean N^^'^A(t{j,m)) 
for any individual sample. Thus, the sum of all samples in the 
jth bit interval, iV^ 6 ;i^[j] = Em=i (^(j. "*))’ i® also 

a Poisson RV whose mean is the sum of the means of the 
individual samples, i.e., = Em=i w)). 

Due to the independent movement of the molecules, node 
q observes molecules that were emitted at the start of the 
current or any prior bit interval. As a result, the number of 
molecules observed within Vg in the jth bit interval due to the 
transmission of sequence is also a Poisson 

RV with mean 

J M 

b1 = E W E ((7 - ^)T + tm). (6) 

2=1 m—1 

where NaII] is the number of molecules released at node n 


where k[j] is the amplification factor of node R in the jth bit 
interval. 

At the end of the 7 th bit interval, node D makes a decision 

(SR) 

on its received signal, and node R receives signal 
In the second phase, node S transmits information bit Ws [j + 
1], and relay node R emits A^Aab + 1] = ^b + l]Ef’A?b1 
molecules, which convey the information regarding Ws[j]. At 
the end of the (j -f l)th bit interval, node D makes a decision 
on Ws[j], WdIj + !]■ Thus, node D receives the jth bit with 
a one bit interval delay. The total duration of transmission for 
a sequence of length L is [L -\- 1)T. 

B. Expected Error Probability 

Given kFs[j], an error occurs in the (j -f l)th bit interval 
if WdIj + 1] ^ Wjsb]- Thus, the error probability of the jth 
bit, F’ebWsb]] written as 

Pe [j I Ws [j]] = W{Ws b] ^ Wn b + 1]) • (9) 

Let us assume that is given. Then, the error proba¬ 

bility of the jth bit when Wsb] = 1 Wjsb] — 0 can be 
written as 


at the beginning of the ith bit interval. 

Pe[j\Ws[j] = = 

Pr(ESb + 1] < UWs[j] = 1, wE), 


The cumulative distribution function (CDE) of the weighted 
sum in the jth bit interval is given by |[T3] Eq. (38)] 

( 10 ) 

pr(Er;A b] <C9iw^„.i) =exp(-EZAb]) 

and 


(YUiil)" 

-< E ^ E (7) 

Pe[j\Ws[j] = 0,W^S~l] = 

pKEE'^^b +1 ] > ^d\ws[j] = 0 , wE), 

( 11 ) 


u ;=0 


III. PERFORMANCE ANALYSIS OE AE RELAYING respectively. Hence, the expected error probability is given by 


In this section, we evaluate the expected error probability 
of the two AE relaying protocols, i.e., variable-gain and fixed- 
gain AE relaying, respectively. 


i-ii 


Peim^Sp] = PlPe[j\Ws[j] = l,W^s~ 

+ PoPe[j\Ws[j]=0,W^S~l] 


( 12 ) 






In 1^, for given [i], i e {1, • • • , j}, [j + 1] 

is a Poisson RV with mean 

J M 

m] =E w E Pob:Miu-i)T+t^)^ 

2=1 m—1 

(13) 

where, in turn, Ef W ^ Poisson RV. Thus, the 

conditional probability in ( [T0| ) can be evaluated as 

Pe{j\Ws[j] = 1 ,WE] = 

oo oo 

i: - ■ i: pr (<£>1,+ii <foiivif «iii=7., ■ ■ ■. 

71=0 7j=0 

Eff b1 = 7„ l^sb1 = 1, wE) X Pr (ECf [1] = 71, 


• • • ,Emi b1 = 7,Ksb1 = 1, , (14) 

where ji, i G {I,- - • ,j}, is one realization of Poisson RV 


N. 


(S.fl)r 


Ob, A, L*] for given {Ws[j\ = 1, W^.;^ }. The first probability 
inside the sums in d is the conditional CDF of Poisson 
RV + 1] which can be evaluated based on 1^ 

given N\fj'A 2 b + !]• Thi^ time-varying mean, in turn, can 


be evaluated based on (|13|l after substituting Ef El*] 'vi^h 


fc[* + l] 7 i. The second proBability inside the sums in (14 1 is the 
conditional joint probability mass function (PMF) of Poisson 
E^Ai [*],* G {1, • • • , j}, for given {Ws[j] = 1, W^}, 
where due to the independent observation of molecules in each 
bit interval, we have 


Pr 


(EfE[i]=7i,--- b]=7,E5bi = i,wE 


n 


exp(-iV^f;f j [i]) 


(15) 


where N\^’ [i] can be evaluated based on ( 61 after substitut¬ 
ing A^a[*], n, and q with NA-iWs\i\, S, ana i?, respectively. 
Analogously, the conditional probability in d can be eval¬ 
uated by considering the complementary probability of ( [l4| ) 
after substituting lTs[j] = 1 with tTs[j] = 0. In ( [l4] l, there 
are infinitely many realizations for each Poisson RV E [*]. 
In order to keep the complexity of evaluation low, we consider 
only one random realization of each Poisson RV K 


with mean N 


Af(S,R), 


(S,fl)r 

ob,Aj_[ 


ob,Ai [ 


for a given for evaluation of (14 1 . 


Our simulation results in Section |lVj confirm the accuracy of 
this approximation. The expected average error probability, Pg, 
is obtained by averaging ( [T2| over all possible realizations of 
and all bit intervals. 


C. Optimal Amplification Factor 

In order to maximize the performance of variable-gain AF 
relaying, we derive a closed-form expression for the optimal 
amplification factor, kopt[-], at the relay node for minimization 
of an approximation of the expected error probability of this 
network. Specifically, we approximate the random observation 
of molecules at node R by the average number of observed 




molecules for evaluation of (m|, i.e., in |l3| . 
substituted with its mean a[ [z]. This approximation allows 
us to approximate [j +1 ] in 0 by a Poisson RV wi 

a deterministic mean given as 


with 


j M 

ES[j+i] w Y.Pof:A2 i{j-^)T+tm). 

i—l m—1 

(16) 

We now derive an expression for kopt[‘], for minimization of an 
approximation of (12i. To this end, by substituting 
from ( 61 , [j + 1] can be expressed as 


ESb'+i]E E 

2=1 LJ—l 

M M 

X Y.pSSii^-^)T+tl) (17) 

1 — 1 m—1 


Eq. (17 1 can be rearranged to express -f 1] as a 

functionof the most recent bit, lFs[j], and all prior transmitted 
bits {Ws[i],iG {I,-- - ,7-1}), as 


ESb'+i]= iva, E ^[*+1] E 

2=1 LO—1 

M M 

X EES E pSM{u-i)T^m)+kmNA, 

1—1 m—1 

j-1 M M 

X E w^sMEES((j-^)r+ib EESb™) 

uj—l 1—1 m—1 

M M 

+ k[j + i]iv^,iE 5 b]EESbb ( 18 ) 

1=1 771=1 

where the first term on the right hand side represents the ISI 
produced in the second hop, and the second term represents 
the expected number of molecules observed in Vd due to the 
amplification of molecules originating from previous bit inter¬ 
vals and observed at node R in bit interval j, Ws[i]:i < j (we 
refer to this sum as the amplified ISI of the first hop). The third 
term is the expected number of molecules observed due to the 
transmission of the most recent bit W 5 [j]. Let us assume that 
the fc[z -f 1] are given. Considering ( [Tsj l and the use of a fixed 
threshold at node D, we observe that decreasing k[j -f 1] 
reduces the overall effect of ISI, and increases the probability 
of miss detection, i.e., ^r{Wo{j + 1] f W^sbllH^sb] = !)■ 
On the other hand, increasing k[j + 1] enhances the ISI 
which in turn increases the probability of false alarm, i.e., 
Pr(lV£)b + 1] 7^ ff"sb]|l^sb] = 0)- Thus, the expected error 
probability in ( p^ can be minimized by optimizing k[j + 1]. 
The optimal k[j + 1] is provided in the following proposition. 

Proposition 1: Given W^Tf, the optimal amplification factor 
at relay node R at the beginning of the (j -1- l)th bit interval. 













kopt[j + 1]^ which minimizes the approximate expected error 
probability of the jth bit, can be approximated as 


TABLE I 

SYSTEM PARAMETERS USED IN SIMULATIONS 


^opt [j 1 ] — 


'"Y 


PiBiU 


Po Bo U 


Y^ivi p( 0 ,u)/, p{R,D),, 

1^1=1 ^ob,Ai 


(19) 


where [•] is the nearest integer, In(-) is the natural logarithm, 
and Bb[j + 1], 6 S {0,1} is obtained using (17l as 


Parameter 

Symbol 

Value 

Probability ot binary 1 

Pi 

03 

Receiver distance 

Xd 

500 nm 

Length ot transmitter sequence 

L 

M) 

Radius of relay node H 

tr 

45 nm 

Radius of node D 

ro 

45 nm 

Diffusion coefficient |14|, |i3] 

Dxf 

4.365 X lO-iOsX 

S 


j-i 


M 


’ LU^l 1^1 

M M 

X Pobfjitrn) + Ws[j] ^20) 


Z =1 


m—1 


Proof: N^A;S [j + 1] is approximated as a Poisson RV, 
and hence has a discrete CDF which makes the optimization 
of k[j + 1] difficult. Fortunately, this CDF can be well 
approximated by a continous regularized incomplete Gamma 
function as 0 Eq. (1)] 






mn) 


( 21 ) 


where r(s,p) is the incomplete Gamma function given by 
r(s,p) = e~*P~^dt 117 Eq. (6.5.3)]. The Gamma 
function, r(s), is a special case of the incomplete Gamma 
function with p = 0. Using this approximation, we take the 
partial derivative of (12i with respect to k[j + 1], solve the 
resulting equation for k\j + 1], and round the result to the 
nearest integer value kopt[j + 1], which yields ( [T^ . ■ 

In the remainder of this subsection, we differentiate between 
two types of variable-gain AF relaying. 

Type 1: In this case, we assume that relay node R adjusts its 
amplification factor in the jth bit interval to kopt[j]- However, 
as shown in (19 1 , the evaluation of kopt[j] requires knowledge 
of the bits transmitted by node S. In order to cope with this 
uncertainty, relay node R compares its received message in 
the jth bit interval, a decision threshold 

to estimate VUs[j]. We denote this estimate by Wjj[j], and for 
compactness the estimated information bits up to the current 

^ j — \ 

bit are denoted by . 

* i-i 

Thus, given relay node R estimates the optimal 


amplification factor in the jth bit interval, kopt[j + 1], via (191 


after substituting with , and releases [j + 1] = 


J-i 


(S R) 

opt[j + 1] X ^2 molecules at the beginning of the 


k, 

next (j + l)th bit interval. 

Remark 1: The complexity associated with Type 1 variable- 
gain AF relaying is the same as that of DF relaying, since for 
estimation of kFs[j], i.e., TUrY], relay node R has to decode 
the received signal in each bit interval, which is not desirable. 

Type 2: To reduce the computational complexity of Type 1 
variable-gain AF relaying, we substitute kept [j + 1 ] by kopt [j + 
1 ], where kopt[j + 1] is obtained by averaging (19i over all 


possible realizations of the sequence i.e.. 


kopt [j + 1] — 


1 


W 


^ ^ ^opt [j + 1] I 


( 22 ) 


w 


where IW is a set containing all possible realizations of 
and I yV I denotes the cardinality of set W. Hence, for Type 
2 variable-gain AF relaying the amplification factor in each 
bit interval can be computed offline. In Section IV we show 
that kopt\'] converges to an asymptotic value after several bit 
intervals. 

The overall end-to-end expected error probability of both 
Type 1 and Type 2 variable-gain AF relaying can be evaluated 
via (12 1 after substituting k[i+l\,i G {1, • • • , j}, with kopt[i+ 
1 ] and kopt\i + 1], respectively. 


D. Fixed-Gain AF Relaying 

In fixed-gain AF relaying, the relay node R adopts a 
bit-interval-independent, fixed amplification factor k for for¬ 
warding the signal received from node S. Specifically, the 
amplification factor is denoted by kopt, where kopt is obtained 


by averaging (19 1 over all possible realizations of W' 
over all bit intervals, i.e.. 


j-i 

S;1 


and 


1 , 

~ r \ \ ’ ^23) 

where set Af contains all possible realizations of the sequence 

Fixed-gain AF relaying has a lower real-time computational 
complexity in comparison to both proposed variable-gain AF 
relay protocols, since no bit detection and no adjustment of 
the amplification factor {kopt[']) are required at the relay node. 
Also, no prior knowledge of time-varying amplification factors 
in each bit interval (fcopt[’]) is required. The overall end-to-end 
expected error probability of fixed-gain AF relaying can be 
obtained from ( [I^ after substituting k[i -f l],z G {1, • • • , j}, 
with kopt- 


IV. NUMERICAL RESULTS 

In this section, we present simulation and analytical results 
for evaluation of the performance of the proposed relaying 
protocols. We adopted the particle-based stochastic simulator 
introduced in IB- In our simulations, time is advanced in 
discrete steps of tg, i.e., the time between two consecutive 
samples, where in each time step molecules undergo random 
motion. The environment parameters are listed in Table |I] 































Fig. 2. Average error probability of 10th bit, given one randomly-chosen 
9-bit sequence, vs. amplification factor k. The first 9 bits are “101101001”. 


In order to focus on a comparison of the performance of the 
different relaying protocols, we keep the physical parameters 
of the relay and the receiver constant throughout this section. 
The parameters that we vary are the decision threshold 
the modulation bit interval T, the amplihcation factor k, and 
the frequency of sampling 1 /fo- 

In the following, we refer to the case when no relay is 
deployed between node S and node D as the baseline case. 
We adopt — 2500 for the two-hop network. For a 

fair comparison between the two-hop case and the baseline 
case, we assume for the baseline case that for transmission 
of information bit “1”, = 2500 -I- 2Na2 molecules are 

released by S, where Na 2 is the average number of molecules 
released by node R for transmission of one bit in the two-hop 
case. 

In Fig.|^ the average error probability of the 10th bit, given 
one randomly-chosen 9-bit sequence, is found as a function 
of amplihcation factor k, for system parameters M = 10, 
T = 400 ^s, and = 20. The simulations are averaged over 
10® independent transmissions of the chosen sequence. We see 
that increasing k increases and decreases the probabilities of 
false-alarm and miss-detection, respectively, which conhrms 
the existence of an optimal k for the minimization of ( 12 i. 
The expected error probabilities are evaluated using ( [TOl l and 
( 111 , when VFsflO] = 1 and Ws[ 10 ] = 0 , respectively, 
and considering two approximations for i.e., Eqs. 

and Although using the approximation leads to an 
underestimation of the expected error probability in Fig. 1 ^ the 
difference between kopt[-] obtained via ( [I^ and kopt[-] found 
via simulation is small (less than 9%). 

In Fig. we show the amplihcation factor for Type 2 
variable-gain AF relaying, kopt[j], as a function of the bit 
interval j, for three different sets of system parameters. 
The results are averaged over 10® independent realizations 
of Wg.]^. For the considered scenario, kopt[j] first increases 
and then decreases to an asymptotic value after several bit 
intervals. This behavior is mainly due to the fact that, after 
several bit intervals, the ISI approaches an asymptotic value; 
see GD- The slight difference between the amplihcation factor 
for hxed-gain AF relaying, kopt, and the asymptotic values of 
kopt[j] is due to the fact that kopt[j] is averaged over all bit 



Fig. 3. Time-varying kopt\j\ as a function of bit interval. 



Fig. 4. Average error probability of a two-hop link as a function of the 
amplification factor k. 


intervals to obtain kopt- 

In what follows, the simulated error probability is averaged 
over 3x10“^ random realizations of the sequence and the 

expected error probability was evaluated via (12 1 , after taking 
into account the modihcations required for each protocol. 
Furthermore, for a fair comparison of the performance of 
the network for different bit intervals, we assume that the 
frequency of sampling, to, and the number of samples per 
bit interval, M, are both independent of T, i.e., for any T 
the samples are taken at times fm = {20,40,60,..., 200} ps 
within the current bit interval. 

In Fig. we evaluate the average error probability of a two- 
hop network with hxed-gain AF relaying to assess the accuracy 
of the optimal value kopt- We evaluated for different 
system parameters, i.e., M = {10,20}, T = {400, 600} /is, 
and = {10,20}. We observe that, by doubling 
and keeping the other parameters constant, the optimal k is 
approximately doubled which is in agreement with ( |T^ , when 
Pi — Pq = 0.5. We can also see that increasing M and T 
decreases and increases the optimal k, respectively. Finally, 
we note the excellent match between the optimal k observed 
via simulation and that derived analytically. 

Fig. 0 shows the average error probability of hxed-gain 
AF relaying and the baseline case as a function of detection 
threshold for two different sets of system parameters, i.e., 
{M = 10,T = 400ps,fc = 200} and {M = 10,T = 
600 fis,k = 250}. The results show that hxed-gain AF 
relaying improves the overall performance of the network. We 

































V. CONCLUSION 



Fig. 5. Average error probability of fixed-gain AF relaying and the baseline 
scheme vs. the detection threshold 


In this paper, we considered a three-node network where 
a nano-relay is deployed between a nano-transmitter and a 
nano-receiver. We assumed that the nano-transmitter emits 
multiple random bits, and proposed two relaying protocols, 
namely fixed- and variable-gain AF relaying. Furthermore, we 
derived closed-form expressions for the optimal amplification 
factors at the relay node for minimization of the expected error 
probability of the network. We showed via simulation and 
analysis that AF relaying improves the overall performance 
of the network. In particular, fixed-gain and Type 2 variable- 
gain AF relaying are attractive for implementation when the 
relay node has limited computational capabilities. 
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Fig. 6. Average error probability of fixed- and variable-gain AF relaying, 
DF relaying, and the baseline case as a function of modulation bit interval T. 


also see that increasing T improves the performance of fixed- 
gain AF relaying, since increasing T reduces the effect of ISI 
which, in turn, decreases the effect of amplified ISI at the 
destination node. Furthermore, a comparison of the results in 
Fig. I^and Fig. [^reveals that the optimization of the detection 
thresholds at node D for a given k is equivalent to optimizing 
k for a given 

In Fig. 1^ the average error probability of a two-hop network 
for fixed- and variable-gain AF relaying, the DF relaying 
scheme from o), m), and the baseline case are evaluated as 
a function of modulation bit interval T. kopt, kopt [■], and kopt [•] 
are applied for fixed-gain. Type 1 variable-gain, and Type 
2 variable-gain relaying, respectively. The results show that 
all considered AF relaying protocols outperform the baseline 
scheme where the performance gap increases with T. Both 
types of variable-gain AF relaying perform slightly better 
than fixed-gain AF relaying, since the amplification factor is 
adjusted in each bit interval according to the current optimal 
amplification factor based on the expected ISI. Furthermore, 
we observe that DF relaying outperforms AF relaying. This is 
because in DF relaying, the ISI of the first hop is not amplified 
by the relay node. However, the improved performance of DF 
relaying comes at the expense of an increased complexity as 
the relay has to decode the message received from the source, 
which is not necessary for fixed-gain and Type 2 variable-gain 
AF relaying. Hence, these schemes are more suitable for nodes 
with limited processing capability. 
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